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Abstract 

We review some results of applying the higher covariant derivative regularization 
to the investigation of quantum corrections structure in N = 1 supersymmetric 
theories. In particular, we demonstrate that all integrals, defining the Gell-Mann- 
Low function in supersymmetric theories, are integrals of total derivatives. As a 
consequence, there is an identity for Green functions, which does not follow from 
any known symmetry of the theory, in N = 1 supersymmetric theories. We also 
discuss how to derive the exact /3-function by methods of the perturbation theory. 
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1 Introduction. 



Supersymmetry is certainly one of the most prominent achievements of the high en- 
ergy physics. Soon after its discovery it was found that in supersymmetric theories the 
ultraviolet behavior was essentially improved due to some non-renormalization theorems. 
For example, there are no divergences in the N = 4 supersymmetric Yang-Mills theory, 
and in theories with N = 2 supersymmetry divergences are present only in the one-loop 
approximation. That is why such models are very attractive from the theoretical point of 
view. But, possibly, the most interesting fact is an indirect experimental proof of super- 
symmetry existence in the Standard model. It was obtained by precise measuring three 
coupling constants and investigating their evolution by the renormgroup equations. The 
result is that only in the supersymmetric version of the Standard model the coupling con- 
stants coincide at some energy scale, which should follow from Grand Unification theories. 
So, there is no doubt that investigating supersymmetric theories and, in particular, their 
quantum properties, is very interesting. Dynamics of supersymmetric theories is highly 
nontrivial. It is worth mentioning summation of instanton corrections in the N = 2 su- 
persymmetric Yang-Mills theory pQ (see also [21 [3l HJ [5]) or finding a relation between 
the N = 4 supersymmetric Yang-Mills theory and a string theory, compactified on the 
manifold AdS$ x S5, known as the AdS/CFT-correspondence [6]. 

In this paper we will consider theories with unextended supersymmetry. Such mod- 
els are especially interesting, because the physics seems to be N = 1 supersymmetric 
at energies of the order 10 3 GeV. Dynamics of N — 1 supersymmetric theories also has 
some interesting features. For example, as a result of investigating instanton contribu- 
tions structure, in Ref. [7] form of the /3-function was suggested exactly to all orders of 
the perturbation theory. This /3-function, called the exact Novikov-Shifman-Vainshtein- 
Zakharov (NSVZ) /^-function is 



a 2 



3C 2 -C(i*)(l- 7 (a)) 



2tt(1 - C 2 a/2n) ' v ' 

where 7(a) is the anomalous dimension of the matter superfield in a representation R, 
and C(R) is defined by 
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tr (T a T b ) = C(R) 5' 



(2) 



Such a /3-function has not yet been derived by methods of the perturbation theory. Its 
numerous verifications by explicit calculations up to the four-loop approximation were 
made in Refs. jH [9j [10] • The authors calculated the /3-function, defined by divergence 
in the MS'-scheme. The main result is the following: if a subtraction scheme is tuned 
by a special way, then it is possible to obtain the exact NSVZ /3-function. Nevertheless, 
there is an open question, in which scheme one can obtain this /3-function. The answer to 
this question is, in particular, given in this paper. Moreover, we show that there are new 
identities, relating Green functions in some N = 1 supersymmetric theories, which do not 
follow from known symmetries of the theory. Possibly, this assumes existence of some 
new invariances. Note that possibility of their existence was discussed, for example, in 
finite N = 1 supersymmetric theories [11]. Thus, the dynamics of iV = 1 supersymmetric 
theories is highly nontrivial and deserves further investigation. An important constituent 
of this investigation is a regularization [12]. The matter is that the dimensional regular- 
ization [13] breaks the supersymmetry and is not convenient for studying supersymmet- 
ric theories. Most calculations were made with the dimensional reduction [14|. It is a 
modification of the dimensional regularization, which does not break the supersymmetry 
explicitly. But this regularization appears to be inconsistent [15], and its using can lead 
to some artifacts. A consistent regularization, which does not break the supersymmetry, 
is the higher covariant derivative regularization [16] . Despite of these attractive features, 
its using is technically complicated. That is why before recently this regularization was 
applied only once, for the one-loop calculation in the (non-supersymmetric) Yang-Mills 
theory [T7] . Taking into account comments, made in subsequent papers [18], [191 120] , the 
result of the calculation coincided with the the standard expression for the one-loop j3- 
function (although in original paper [TT] the authors affirm that it is not so). (For other 
applications of higher derivatives see, for example, [21] and the references therein.) 

In this paper we will try to demonstrate that using the higher covariant derivative reg- 
ularization allows revealing some interesting features of the quantum correction structure 
in N = 1 supersymmetric theories, which was not known earlier. 

This paper is organized as follows. 

In Sec. [2] we recall basic information about the N = 1 supersymmetric Yang-Mills 
theory, the background field method, and the higher covariant derivative regularization 
in supersymmetric theories. Then, in Sec. [3] we present results of explicit calculations, 
made in supersymmetric theories with the higher covariant derivative regularization, and 
analyze their features. After that, in Sec. H]we try to explain them using the Schwinger- 
Dyson equations and Slavnov-Taylor identities. It turns out that in order to explain 
results of explicit calculation and to obtain the exact NSVZ /3-function it is necessary 
to propose the existence of a new identity for Green functions. Different forms of this 
identity are discussed in Sec. [5j In order to verify this identity, some calculations are 
made in the three- and four-loop approximations. They are described in Sec. [6j Structure 
of quantum corrections in supersymmetric theories, obtained with the higher derivative 
regularization, is discussed in the Conclusion. 
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2 N = 1 supersymmetric Yang— Mills theory, back- 
ground field method, and higher derivative regu- 
larization 



In this paper we consider the N = 1 supersymmetric Yang-Mills theory with matter 
fields, which is described in the superspace by the action 

S = ^Retr J d 4 xd 2 9W a C ab W b + ^ J d 4 xd 4 6 (Ve 2y + ^ + e~ 2Vt pj + 
+-m [ d 4 xd 2 9ft0 + -m [ d A xd 2 0^ + <p*. (3) 



2 7 2 

Here and are chiral matter superfields and V is a real scalar superfield, which contains 
the gauge field as a component. The superfield W a is a supersymmetric analogue of 
the gauge field stress tensor. It is defined by 



w a = Id 2 



e- 2V D a e 2V 



(4) 



where D a and D a are the right and left supersymmetric covariant derivatives respectively. 
In our notation the gauge superfield V is decomposed with respect to the generators of 
a gauge group T a as V = e V a T a , where e is a coupling constant. The generators of the 
fundamental representation we will denote by t a . They are normalized by the condition 

tr(t a t b ) = ^6 ab . (5) 
Action (jHJ) is invariant under the gauge transformations 

-> e* A 0; $ - e~ iAt 0; e 2F - e iA+ e 2l V 4A , (6) 

where A is an arbitrary chiral superfield. Such a transformation law means that if the field 
is in the representation R of the gauge group G, then the field <fi is in the representation 
R, conjugated to R. 

Note that theory is not the most general renormalizable supersymmetric model. 
In principle, it is possible to consider a case, in which matter superfields are in an arbi- 
trary representation of a gauge group (instead of R + R), and add terms, cubic in the 
matter superfields. However, calculations in model are simpler. That is why here we 
investigate only such a theory. 

For quantization of this model it is convenient to use the background field method. 
The matter is that the background field method allows calculating the effective action 
without manifest breaking of the gauge invariance. In the supersymmetric case it can 
formulated as follows [221 [23] : Let us make a substitution 



e 2v ^ £ 2v> = e ^ e 2v e n (7) 

in action where f2 is a background scalar superfield. An expression for V is a 
complicated nonlinear function of V, CI, and Cl + . We do not interested in explicit form 
of this function: 
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v' = v'[v,n]. (8) 

(For brevity of notation we do not explicitly write the dependence on f2 + here and below.) 
The obtained theory will be invariant under the background gauge transformations 



(j) -> e iA (j); (j) -> e~ lA 0; V — > e lK Ve- lK ; 

e « _> e ^e"e- jA ; e" + - e iA W^, (9) 

where K is an arbitrary real superfield. However, there is one more invariance. In order 
to construct it, we define first the background chiral covariant derivatives 

D a = e- n+ D a e n+ - D a = e n D a e- n . (10) 

Acting on some field X, which is transformed as X —>■ e iK X, these covariant derivatives 
are transformed in the same way. It is also possible to define a covariant derivative with 
the Lorentz index 



D^- l -{Cr) ah {D a) D b ] ) (11) 



which will have the same property. It is easy to see that the action is also invariant under 
the quantum transformations 



e «0 ^ e ^ e n . e -n^ ^ e -* e - nt 0; e 2V - e lA+ e 2V V a ; 

o, -»• n, n + -> n+ (12) 

where A is an arbitrary background chiral field, which satisfies the condition D\ = 0. 
Such a superfield can be presented as A = e n Ae _fi , where A is a usual chiral superfield. 
It is easy to see that after substitution <£7§ action (j3J) will be 

S = _L t rRe J d A x d 2 9W a W a - -^ti Re J d A x d A 9 lQ>{e~ 2V D a e 2V ^jW a + 



+ (e- 2V D a e 2V )D 2 e- 2V D a e 2V 



xe-"^) d^xdHft^+^m / d 4 x d 2 6 4 >*, (13) 

where 



W a = | e »£H (e^e"" D a e" e»j e~". (14) 

Action of the covariant derivatives on the field V in the adjoint representation is defined 
by the standard way. 

It is convenient to choose a regularization and gauge fixing so that invariance (Q will 
be unbroken. First, we fix a gauge by adding 



5 



s, 



gf 



32e 



^tr J d A xd A 6 (vD 2 D 2 V + VD 2 D 2 V^j 



(15) 



to the action. In this case terms quadratic in the superfield V will have the simplest form: 

^trRe J d A xd A eVD\V. (16) 
The corresponding action for the Faddeev-Popov ghosts S c is written as 



S c = itr j d A xd A 9 !^(c + c + )V (c + c+) + cth^(c - c+) | 



(17) 



The superfield V in this expression is decomposed with respect to the generators of the 
adjoint representation of a gauge group, and the fields c and c are the anticommuting 
background chiral fields. 

Moreover [22], the quantization procedure also requires adding the action for the 
Nielsen-Kallosh ghosts 



S B = ^Ltr J d A xd A 6B + e n+ e n B, (18) 

where B is an anticommuting chiral superfield, and the background field should be de- 
composed with respect to the generators of the adjoint representation of a gauge group. 
Because the fields B and B + do not interact with the quantum gauge field, they contribute 
only to the one-loop (including subtraction) diagrams. It is important to note that the 
factor 1/e 2 in action ( JTBl is the same as in action for the gauge fixing terms (1151) . 

The gauge fixing breaks the invariance of the action under quantum gauge trans- 
formations (j!2p . but there is a remaining invariance under BRST-transformations. The 
BRST-invariance leads to Slavnov-Taylor identities, which relate vertex functions of the 
quantum gauge field and ghosts. 

For regularization we use the following method: Let us add the term 



l r (D 2 ) n+1 

Sa = ^trRe J d A xd A 9V K -^V (19) 

to action (|T3|) . Then the invariance under both supersymmetry transformations and trans- 
formations (Q is unbroken. Therefore, the effective action, calculated with the background 
field method, is invariant under both supersymmetry and background gauge transforma- 
tions. The described way of regularization is a bit different from the method, proposed 
in Ref. |24j . The difference is in the form of the term, which contains higher covariant 
derivatives. In the method, considered here, it breaks the BRST-invariance, but the form 
of terms, quadratic in the quantum superfield V, is simpler. This simplifies calculations 
in a certain degree, while all particular features of the higher derivative regularization 
are the same in the both cases. However, because the higher derivative term breaks the 
BRST-invariance, it is necessary to use a special subtraction scheme, which cancels non- 
invariant terms and guarantees fulfilling the Slavnov-Taylor identities in each order of 
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the perturbation theory. Such a scheme was proposed in Refs. [25, 26J and generalized 
to the supersymmetric case in Refs. [27J EE] • With the background field method such a 
scheme is simpler, because the background gauge invariance guarantees, for example, the 
transversality of the two-point Green function of the gauge field. Nevertheless, additional 
subtractions should be made for Green functions, containing the ghost fields. 
Let us construct the generating functional as follows: 

Z[J, f2] = J Dfi exp f^iS + iS A + iS gf + iS gh + iS s + iS<p + 

+i Jd 4 xd 4 e (j + j[n]j(y'[v,n] - v)}, (20) 

where S g f is gauge fixing terms (fl5T) and S g h = S c + Sb is the corresponding action for 
the Faddeev-Popov and Nielsen-Kallosh ghosts. In all expressions the coupling constant 
e should be substituted by the bare coupling constant 6q. Ss denotes terms with sources 
for chiral superfields. In the extended form they are written as 

S s = J d 4 xd 2 e(j t (f> + ] t ^j + J d 4 xd 2 9(j + ^+j + 4>*y (21) 
Moreover, in generating functional (|20|) we introduce the additional sources 

S H = ~J d A xd 4 6 (0 o + e" + e 2V 'e n + ^ e n+ e 2V e n ^ + 

+^e~ n *e- 2Vt e- nt 4> + 4> + e- n *e- 2Vt e- nt 4>o), (22) 



where (j) , 0q , (fi and (pQ are arbitrary scalar superfields. In principle, it is not necessary 
to introduce the term in the generating functional, but the presence of the parameters 
4>q is highly desirable for investigating the Schwinger- Dyson equations. The superfield V 
is defined by 

e 2V = e n+ e n (23) 

and J[fi] is a so far undefined functional. A reason of its introducing will be clear later. 
The functional integration measure is written as 

Dfi = DV DcDcDB D<j>D<j>. (24) 

In order to understand how generating functional fl20l is related with the ordinary 
effective action, we perform the substitution V — > V. Then we obtain 

z\j,n] = exp | -i J d 4 xd 4 e (j + j[n])v}z J + J\n],n , (25) 

where 

UJ, 11] = /^ exp { tS + ,S A + iS„ + tS , h + lS$ + ,5,, + i J A d'e JV\ (26) 
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If the dependence of S, S\, S g j, S g h, and Sfa on the arguments V, f2, and f2 + were 
factorized into the dependence on the variable V, Z would not depend on ft and S7+ 
and would coincide with the ordinary generating functional. This really takes place for 
action ([3]). However, in the term with the higher derivatives, in the gauge fixing terms, 
and in the ghost Lagrangian such factorization does not occur. Therefore, Zq actually 
differs from the usual generating functional. 

Using the functional Z[J, fl,j) it is possible to construct the generating functional for 
the connected Green functions 



w[j,n,j] = -i\nz[j,nj] = 

= - [ d 4 xd 4 9(j + j[n])v + w 



and the corresponding effective action 



The sources should be expressed in terms of fields using the equations 



(27) 



r[v, n, 0] = - J d 4 xd 4 e (jv + j[o]v) + w j + J[Q], n, j - 
- [ d 4 x d 4 e jv- [ d 4 x d 2 e (j* 4>+j t $)- [ d 4 x d 2 9 (j + </>* + j V ) . (2* 



v = jjW[j, o, j] = -v + jjW j + n, 3 



^w[j,n,j] = -^w 
03 03 



e.t.c. 



Substituting these expressions into Eq. (128]) . we write the effective action as 



(29) 



r[^,o,, 



d 4 xd 4 e ( j[n]v + Jyi w o 3, J + JM,nj 



SJ 



d 4 xd 2 6<j)—W 



j+j[n],n,j 



— ( similar terms with 



!>, and 



(30) 



Let us now set V = 0, so that 



(31) 



We also take into account that the invariance under background gauge transforma- 
tions (Q essentially restricts the form of the effective action. If the quantum field V in 
the effective action is set to 0, then the superfield K will be present only in the gauge 
transformation law of the fields Q and Q + , the only invariant combination being expres- 
sion (I23p . (It is invariant in a sense that the corresponding transformation law does not 
contain the superfield K.) This means that in the final expression for the effective action 
we can set 
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v. 



(32) 



In this case the effective action is 



r[o,v,0] = w / o J + J[v],v,j - d 4 xd 4 e(j + j[v 



- I d A xd 2 e <p^-w 
03 



J + J[V),V,j 



— ( similar terms with 



J + J[V],V,j 
and c6 



(33) 



Note that this expression does not depend on form of the functional J[fl]. In particular, 
it can be chosen to cancel terms linear in the field V in Eq. (120]) . Such a choice will be 
very convenient below. 

If the gauge fixing terms and the terms with higher derivatives depended only on 
V, expression (I3"3"|) would coincide with the ordinary effective action. However, as we 
already mentioned above, the dependence on V, ft, and fl + is not factorized into the 
dependence on V with the proposed method of regularization and gauge fixing. According 
to Ref. [291 [30], the invariant charge (and, therefore, the Gell- Mann-Low function) is 
gauge independent, and the dependence of the effective action on gauge can be eliminated 
by renormalization of the wave functions of the gauge field, ghosts, and matter fields. 
Therefore, for calculating the Gell-Mann-Low function we may use the background gauge 
described above. We note that if this gauge is used, the renormalization constant of the 
gauge field is 1 due to the invariance of the action under transformations (Q. 

Nevertheless, generating functional (120]) is not yet completely constructed. The mat- 
ter is that adding the term with higher derivatives does not remove divergences from 
one-loop diagrams. To regularize them, it is necessary to insert the Pauli-Villars determi- 
nants in the generating functional [31J. The Pauli-Villars fields should be introduced for 
the quantum gauge field, ghosts (including the Nielsen-Kallosh ghosts) and the matter 
superfields. Constructing them we will at once use condition (|32|) . So, we should insert 
in the generating functional the factors 



(34) 



II(det Py(V,V,M,) 

i 

in which the Pauli-Villars determinants are defined by 



det PV(V, V, M) 



DVpvDcpvDcpvDB PV D(j)pvD(f>pvexip(iSpv), (35) 



where the action for the Pauli-Villars fields is 



S PV = trRe / d?xd?eV PV 



1 WA, + y 2 M v 

e e 



V PV + 



+ -tr / d xd 6(cp V + Cp V )V (c PV + cj, v ) + cth V(c PV - cj, v ) 
4 ( ^M c tr / d A x d 2 9 c PV c PV + h.c.) + tr / d 4 x d A 6 B+ v e 2V B r \ 
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74^4^+ ye V e 2V e V 0py + 



+tr (i / d ^d 2 9M B B 2 PV + h.c}j + \z J d 4 xd 4 9 

+0+ y e- vt e- 2V "e- vt Py ) + Q J d 4 x d 2 9 M^ PV (j> PV + h.c.J , (36) 

and Z is a renormalization constant for the matter superfields. The Grassmanian parity of 
the Pauli-Villars fields is opposite to the Grassmanian parity of usual fields, corresponding 
to them. The coefficients q in Eq. (|34p satisfy the conditions 



5> = 1; £c^ = 0. (37) 

i i 

Below, we assume that Mi = e^A, where dj are some constants. Inserting the Pauli- 
Villars determinants allows cancelling the remaining divergences in all one-loop diagrams, 
including diagrams containing counterterm insertions. (This is guaranteed because the 
masses of the gauge field and Nielsen-Kallosh ghosts are multiplied by the renormalized 
coupling constants, and the other terms are multiplied by the bare ones. This will be 
discussed later in more details.) 

In this paper we will calculate the Gell-Mann-Low function, which is determined by 
dependence of the two-point Green function on the momentum in the limit m — > 0. That 
is why we consider the massless case and write terms in the effective action, corresponding 
to the renormalized two-point Green function, as 

rf = -y^tr j -0^d 4 e\(-p)d 2 n 1/2 V(p)d~\a^/p), (38) 

where a is a renormalized coupling constant. The Gell-Mann-Low function, denoted by 
/3(a), is defined by 

p(d(a,p/p)J = ^jd(a,/i/p). (39) 

The Gell-Mann-Low function is scheme independent. To prove this, we will use the 
following statement as a starting point: If we fix a normalization point fi <C A and 
impose in this point the boundary condition for the renormalized two-point Green func- 
tion d{pj p, = 1), then the two-point Green function is uniquely determined and does 
not depend on both renormalization and regularization. For example, if two different 
regularizations (or renormalization schemes) are used, then 

di(aiOx),^) =d 2 (a 2 M,^), (40) 

where aj(/i) and di are the renormalized coupling constants at the scale \i and the renor- 
malized two-point Green functions, obtained in the first and in the second regularization 
respectively. Setting p = n in Eq. (HOI) , it is possible to find the dependence a.\{a-2). 
Therefore, two different regularizations differ in a finite renormalization of the coupling 
constant. We note that such a renormalization can be gauge dependent and cause the 
gauge dependence of the effective action divergent part. However, the Gell-Mann-Low 
function, which we will calculate below in this paper, does not depend on such a finite 
renormalization, because (setting x = Inp/pt) 



10 



/3i (di(ai, x)J = — di(ai, x) = -^d 2 (a 2 , x) = (3 2 [d 2 (a 2 , x)J = (3 2 (di(a:i, x)J . (41) 

Therefore, the Gell-Mann-Low function is independent of the regularization. In particu- 
lar, a regularization can break the BRST-invariance if the renormalized effective action is 
obtained by subtractions, restoring the Slavnov-Taylor identities. 

The anomalous dimension is defined similarly. First we consider the two-point Green 
function for the matter superfield in the limit m — > 0: 



^ 2) = 4 J J^ d4 ° (^ + (-P^)^ 9 ) + ^(-P^)^ 9 )) ZG ( a ^/P)^ (42) 

where Z denotes the renormalization constant for the matter superfield. Then the anoma- 
lous dimensions is defined by 

7 (rf(a,/i/p)) = -—\nZG(a,v/p). (43) 

3 Calculation of quantum corrections with the higher 
derivative regularization 

3.1 Supersymmetric electrodynamics 



Calculation of quantum corrections with the higher derivative regularization was made 
first for the N — 1 supersymmetric electrodynamics. The matter is that in the Abelian 
case the term with higher derivatives is simpler. Really, the superfield W a is gauge invari- 
ant in the Abelian case. Therefore, the higher derivative term contains usual derivatives, 
instead of covariant derivatives. This essentially simplifies the Feynman rules. For the 
considered model in the Abelian case they can be formulated as follows: 

1. External lines gives the factor 

n / / J^k> ■ • (2"> 4 <x>-)- < 44) 

where the index E numerates external momentums. 

2. Each internal line of the superfield V corresponds to the propagator 

(45) 



(A; 2 + iO)(l + (-l)™P«/A 2ri ) 

3. Each internal line — 0* or — 0* corresponds to the propagator 

1 



D 2 D 2 8\e 1 - 9 2 ). (46) 



16(A; 2 + iO) 
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4. The Pauli-Villars fields are present only in closed loops. Each internal line $ — $* 
or $ — $* corresponds to the propagator 



- 16(*»-M? + iO) 5W( *-* ) - (47) 
and each internal line $ — $ or $* — $* corresponds to 

U and -— V; 2 ■ - DH\6 x -e 2 ) (48) 



4(Jfe 2 - M 2 + iO) v ' 4(P - M 2 + zO) 

respectively. For each loop with the Pauli-Villars it is necessary to add — ^^Q- 

i 

/d 4 k 

6. Each vertex produces integration with respect to the corresponding 9: J d 4 9. 

7. It is necessary to calculate a numerical coefficient for each diagram. 

We first present an expression for the two-loop contribution to the effective action, 
corresponding to the two-point Green function of the matter superfield. It was found 
in Ref. [32j and will be needed later for calculating the three-loop /^-function. This 
contribution can be written in form (l4"2l. where 



G= \ + G 1 + G 2 + 0(a 3 ) (49) 

is a result of calculating the two-loop diagrams without insertions of counterterms on 
matter lines. G\ and G 2 denote the following Euclidean integrals: 



1 = ~ I (2^)4 k 2 {k + p) 2 (l + k 2n /A 2n ) ' ^ ^ 

_ d 4 k d 4 l 4e 4 
" 2 = ' J2^ j2nY k 2 l 2 {k + p) 2 {l + p) 2 (l + k 2n /A 2n ){l + l 2n /A 2n ) + 



+ 



d 4 k d 4 l 4e 4 



(2tt) 4 (2tt) 4 k 2 l 2 (l + p) 2 (k + I + p) 2 (l + k 2n /A 2n )(l + l 2n /A 2n ) 
d 4 k d 4 l Ae 4 {k + I + 2p) 2 

(2tt) 4 (2tt) 4 k 2 {k + p) 2 l 2 {l+p) 2 (k + l+p) 2 {l + k 2n /A 2n ){l + l 2n /A 2n ] 



+ 



r d 4 k je 4 / r d 4 i i 

+ J (2T) 4 k 2 {k + p) 2 {l + k 2n /A 2n ) 2 V 7 (2T) 4 Z 2 (£; + Z) 2 ( ^ 

-^Iwv 4 (i 2 + M 2 m + i) 2 + M 2 ) -M ln -, +h )( 1 + k2n/A2n ) ) • 

Therefore [33], the two-loop renormalization constant for the matter superfield can be 
written as 

Z(e,k/ii) = 1 + - ( In - + Zl 

7T V /i 



y ( In 2 - + Zl In -) - 72 a 2 ( In - + z 2 ) + 0(a 3 ), (52) 
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where 72a 2 is the two-loop anomalous dimension. Here we assume that in order to cancel 
one-loop divergences, the counterterms 



1 / A \ f / d 2n \ 

AS = -_ (in - + bl )Re J dS d 2 9 W a C« b (l + — ) W b + 

+ ( In - + z^j I d 4 x d A 9 (>e 2y + 0*e" 2y 0) (53) 

are added to the action, where bi, z\, and z 2 are arbitrary finite constants. Choosing them, 
we fix a subtraction scheme. It is important to note that ZG is finite by construction up 
to terms of the order a 3 . 

Let us proceed now to the calculation of the Gell-Mann-Low function. Diagrams, 
which define it, were calculated in Ref. |33j. It is necessary to remember that each 
internal line of matter superfields can correspond to both the fields 0, and the Pauli- 
Villars fields. The total three-loop contribution to the effective action can be presented 
in form (138D . where the (sufficiently large) expression for the function d~ x , obtained by 
explicit calculating Feynman diagrams, is presented in Ref. [33]. As a check of the 
calculations we verify cancellation of all noninvariant terms proportional to 

/ (I^^MW-^)- ( 54 ) 
In order to construct the Gell-Mann-Low function we consider the expression 



d 



din A 



d 1 (a ,A/p) 



d 



p=0 



din A 



{A l + A 2 + A 3 + A A ) 



(55) 



p=0 



where A\ is a one-loop result; A2 is a sum of two-loop diagrams, three- loop diagrams with 
two loops of matter superfields, and counterterm diagrams, arising due to the renormal- 
ization of the coupling constant; A 3 is a sum of three- loop diagrams with a single loop of 
matter superfields; and A4 is a sum of diagrams with counterterms insertions on lines of 
matter superfields. Results, found in Ref. [23], can be written as 



dAi 



din A 

dA 2 



din A 



-E 



-2e' 



p=0 



p=0 



-2e< 



d A q 1 d 



(27r) 4 g 2 dq 2 (g 2 + M"2) 2 ' 
d A q 1 d d f d 4 k 



(56) 



(2tt) 4 q 2 dq 2 din A J (2tt) 4 p A + fc2n/A2n 



{k + qy 



Ci 



((k + q) 2 + M 2 )(q 2 + M. 



dH 



(2vr) 4 t 2 (k + t) 2 



dH 



(2vr) 4 (t 2 + M 2 )((k + t) 2 + M 2 



(57) 
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dA, 



dm A 



f d 4 q Id 

p=0 ~ J (27r) 4 ^ dg 2 



d d 



d 4 k dH 



4e 4 



x 



x 



1 



dg 2 din A J (2tt) 4 (2vr) 4 £2 ^ + k 2n / A 2n j /2 ^ + p n / A 2n j 

+E 



X 



2k 2 



1 



1 



2(g + /) 2 (g + A; + /) 2 (g + Z) 2 (g + fc + If 

1 2A; 2 
+ 



(q 2 + M 2 f 



2((g + Z) 2 + M 2 ) (( g + A; + /) 2 + M 2 )((g + /) 2 + M 2 ) 
2M 2 



((g + k + I) 2 + M 2 ) {(q + k) 2 + M 2 )((g + k + I) 2 + M 2 ) 
2 Mr -2 Mf 



(g 2 + M 2 )((g + I) 2 + M 2 ) {{q + I) 2 + M 2 )((g + fc + Z) 2 + M 2 ) 
d 4 q 1 d 



(58) 



p=0 



A 



-72 a In h z 2 

/i 



(2tt) 4 g 2 dg 2 

g 4 



( - ( - ( In - + Zl ) + % ( In 2 - + zi In -) - 

[ \7T\/i / ir z \ fj, [i J 

A \ 2 (6g 4 + 3g 2 M 2 + Mf)M 2 



(g 2 + M 2 ) 2 16tt 4 V // 



( ln u +21 ) 



2(g 2 + M 2 f 
2g 4 



+ 



a A \ d r d 4 k e 2 

7T n Ji +91 ) din A J (2tt) 4 jfc2(i + k 2n /A 2n ) (Jk + g) 2 + M 2 )(g 2 + M 2 ) 2 J' 



(59) 



Thus, we reveal an important feature of the quantum corrections structure: all contri- 
butions to the Gell-Mann-Low function in supersymmetric theories are integrals of total 
derivatives. Really, in the four-dimensional spherical coordinates 

I^i Kf) = ^ Kf = oo) - f{f = 0) )- (60) 

Using this equality it is possible to calculate all integrals, presented above. The result is 



dA 1 
din A 

d 



16tt 2: 



din A 



p=0 

(A 2 + A3) 



d 4 k d 



1 



x 



p=0 8n 2 J (2n) 4 d\nAk 4 (l + k 2n /A 2n ) 2 
d 4 t 1 



1 + k 



2n I A 2 n^ 



X 



E 



(2tt) 4 t 2 (A; + t) 2 
d 4 t 



(2tt) 4 (t 2 + M 2 )((k + t) 2 + M 2 ) 



d 4 k d 4 l d 



4tt 2 J (2tt) 4 (2tt) 4 din A (1 + k 2n / A 2n ) (1 + l 2n /A 2n ) [2k 4 l 4 k 4 l 2 (k + I) 2 



d ( 



167r 2 dlnA 



[-G 1 -G 2 + -G( 



p=0 



16tt 2 



InG. 



(61) 



Because the integrals, defining A 4 , depend only on A/p, taking the limit A — > oo is 
equivalent to taking the limit p — > 0. Hence, taking into account that q = 1, we obtain 
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A A = — V —{ - (-f\n- + zi) (in 2 - + Zl ln-) - 72 a 2 ( In - + zS) + 

p=0 l07T Z { \7Y \ jj / 7Y Z \ jj fj, / \ /J, J ) 

+— (ln- + , 1 ) |+0(a 3 ) = - I ^lnZ + 0(« 3 ). (62) 
Collecting all contributions to the function d" 1 , we obtain 

c£ _1 (e, fi/p) = h In \n(ZG) + finite terms. (63) 

Therefore, there are divergences only in the one-loop approximation. (This was first noted 
in Ref. [34] .) In order to compensate them the bare coupling constant should be presented 
in the form 

-1 = ---\ n --b 1 + 0(a 3 ). (64) 

ao a 7T fl 

Then the final result can be written as 

d-\e, n/p) = - + In - - ln(ZG) + 0(a 3 ) + finite terms. (65) 
a p 

Using this expression it is possible to construct the three-loop Gell-Mann-Low function 
according to Eq. (I3TJ1) . In this case it coincides with the expansion of the exact NSVZ 
/3-function 

2 

P(a) = ^-(l -7(a)) +0(a 5 ). (66) 

It is worth mentioning a difference between the results of calculations, made with 
the higher covariant derivative regularization and with the dimensional reduction. With 
the higher derivative regularization divergences are only in the one-loop approximation, 
while with the dimensional reduction they appear in all loops. According to Ref. [35J the 
difference of the results for divergences arises because with the dimensional reduction a 
contribution of diagrams with counterterm insertions is 0, while with the higher derivative 
regularization it is not 0. In the three-loop approximation the situation is completely 
similar. The results for the sum of diagrams with counterterms insertions differ due to 
the mathematical inconsistency of the dimensional reduction, which was first pointed in 
Ref. [15]. In particular, the straightforward application of the dimensional reduction 
to calculating anomalies gives zero. Due to the same reasons (see Ref. [35]) the sum 
of diagrams with the counterterm insertions, defining the rescaling anomaly, which was 
investigated in Ref. [35] in details, is also with the dimensional reduction: 

( exp (i(Z - 1) i j d 4 x d A 6 [^e 2V 4> + 4>*e~ 2V ^j ) ) = 1. (67) 
With the higher derivative regularization we find 
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16n 



(68) 



This result, obtained with the higher derivative regularization, in a certain degree confirms 
speculations, made in Ref. [37]. According to this paper, the Wilsonian action Sw is 
exhausted at the one-loop, while the effective action T has corrections in all loops. Now 
we see that Sw should be replaced by the usual renormalized action. 

It is important to note that the existence of divergences only in the one-loop approxi- 
mation is not a physical result. The renormalized effective action and the Gell-Mann-Low 
function are the same in both regularizations. The relation between the Gell-Mann-Low 
function and the nonphysical /3-function, defined by divergence, is broken because the way 
of introducing the regularization produces the dependence of the generating functional on 
the normalization point fi at a fixed value cq. 

3.2 Yang— Mills theory without matter fields 

In the previous section we demonstrated that in the N — 1 supersymmetric electro- 
dynamics all integrals, defining the Gell-Mann-Low function in the limit p — > were 
factorized into integrals of total derivatives. Does this fact take place in the non-Abelian 
case? To answer this question in Ref. [38] the two-loop Gell-Mann-Low function was 
calculated for the pure N = 1 supersymmetric Yang-Mills theory (without matter fields) 
with the higher covariant derivative regularization. 

The one-loop /3-function, calculated with the background field method, is well-known 
[22]. Using the higher covariant derivative regularization does not essentially change 
the calculation and its result [20] . Let us mention the typical features. The quantum 
superfield V does not contribute to the one-loop diagrams, because in the corresponding 
diagrams a number of the spinor derivatives D, acting on propagators, is less than 4. 
Really, a result of calculating any two-point diagram is proportional to 



where x and y are the points, to which the external lines are attached. The result is not 
only if the operator P contains 4 spinor derivatives. However, two vertexes can contain no 
more than 2 spinor derivatives, and propagators of the gauge field do not contain spinor 
derivatives at all. Therefore, all one-loop two-point diagrams are automatically 0. The 
one-loop diagrams with the Pauli-Villars fields, corresponding to the gauge field, are 
due to the same reason. Since the higher derivatives do not change a number of spinor 
derivatives in vertexes, the one-loop contribution of the quantum field is also in the 
regularized theory. 

Therefore, the one-loop two-point Green-function of the gauge field is completely de- 
termined by contributions of the Faddeev-Popov and Nielsen-Kallosh ghosts. With the 
regularization and gauge fixing, described above, the ghost Lagrangians do not depend 
on the presence of higher derivative terms. Due to anticommuting, the contributions of 




(69) 



x=y 
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each ghost field have opposite sign in comparison with the contribution of a chiral scalar 
superfield in the adjoint representation of a gauge group. Therefore, in the one-loop 
approximation the Gell-Mann-Low function is 



(5{a) 



3C 2 a 2 
2tt 



(70) 



The effective action in the two-loop approximation is calculated by the standard way. 
It is contributed by diagrams, schematically presented in Fig. [TJ Usual diagrams are 
obtained by attaching to them two external lines of the background gauge field by all 
possible ways. In Fig. [1] a propagator of the quantum field V is denoted by a wavy line, 
a propagator of the Faddeev-Popov ghosts by dashes, and a propagator of the Nielsen- 
Kallosh ghosts by dots. (We note that they contribute only in the one-loop approximation, 
because the Nielsen-Kallosh ghosts interact only with the background field.) 








X- 



Figure 1: Diagrams, contributing to the two-loop /9-function of the N — 1 supersymmetric 
Yang-Mills theory. 

With the higher derivative regularization the propagator of the quantum field is 

q 2 {l + q 2n /A 2n ) ^ 

(in the Euclidean space after the Weak rotation). Feynman rules for vertexes, containing 
two lines of the quantum field V, are also changed. In particular, the vertex with a single 
line of the background superfield V, which has the momentum p, (it is denoted by a bold 
wavy line) is 




-Qk+p^DrfrysDVh. 



(k + p) 



2n+2 



k 2n+r ' 



V A 2n [(k+p) 2 - k 2 



(72) 

and the vertex with two lines of the background superfield V, which have momentums p 
and —p, is 
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A/ 1 

~ (4V9 2 n 1/2 V + p 2 V 2 ) (l + ( n + 1)— ) + — ((2* + p) 2 V«9 2 n 1/2 V+ 



+V^1-L ^ 2 ^2>> Afe+p) 2n+2 -^ +2 (n+l)fc 2 " \ 2 

According to the calculations, the two-loop contribution of the Faddeev-Popov ghosts 
to the Gell-Mann-Low function is that agrees, for example, with Ref. [39] . (Inte- 
grals, defining the two-point Green function, appeared to be some finite constants for the 
ghosts.) 

As we already mentioned, the total two-loop contribution of the two-point diagrams 
to the effective action can be presented in form (I38I) due to the Slavnov-Taylor identity. 
To find the function d^ 1 up to an unessential constant, we differentiate it with respect to 
In A, and then set the external momentum p to 0. Later we will see that the result is a 
finite constant d\\ 



( -d 1 (a,A/p) 



= d v (74) 

p=0 



din A 

Therefore, the function d~ l depends on the momentum logarithmically 

d _1 (a, A/p) = di In — h const. (75) 
P 

Calculating explicitly two-loop diagrams, presented in Fig. [T] (so far without diagrams 
with counterterm insertions), differentiating the result with respect to In A, and, then, 
setting p = 0, we obtain (in the Euclidean space, after the Weak rotation) 



d\ = 8tt ■ 6na 



d 4 k 1 d d 
{2n)*¥dk 2 d\nA 



k) 2 x 



x(l + (q + kY n /A 2n ) 



2(n+ 1)(1 + k 2n /A 2n 



1 - 2n( 



1 + fc 2n /A 2n 



(76) 



It is important to note that taking the limit p — ► is rather nontrivial, because the final 
result can contain infrared divergent terms, proportional to p or p 2 , or terms, proportional 
to p, but giving a finite contribution to d\. However, the calculation shows that all such 
terms are cancelled. Moreover, the sum of diagrams appeared to be a total derivative with 
respect to the module of the loop momentum, so that the integral with respect to d A k, 
which is contained in Eq. ( J76|) . can be easily calculated by Eq. ( l60i) . All substitutions at 
the upper limit are due to the higher derivative regularization, and only the substitution 
at the lower limit is nonzero. Using equations, presented above, we obtain 



d\ = — 6«o" 



d 



^ f q 2 (l+q 2n /A 2n ) 



(77) 



'din A J (2tt) 4 

This integral can be also easily calculated in the four- dimensional spherical coordinates: 
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di = — [ S\z-a 2 4~^ + ? 2 7A 2 T 2 = + g 2 7A 2 T 2 



7T J (2tt)V dq 2K * ' ' 4tt 3 



3a ° (78) 



4tt 2 ' 



(We note that the result does not depend on the regularization parameter n.) Therefore, 
in the two-loop approximation 

«£'(<*, A/p) = - - — In - - ^1 In - + O(a' ), (79) 

where do denotes the function d, calculated without diagrams with counterterms inser- 
tions. 

Therefore, the Gell-Mann-Low function, defined by Eq. (1391) . in the two-loop approx- 
imation is 

^-^-w +0(a4) (80) 

and coincides with the expansion of the exact NSVZ /3-function in the considered order. 
We note that this result does not depend on a possible finite constant in Eq. (!79l) . 

For calculating quantum corrections it is also necessary to take into account diagrams 
with counterterms insertions. Usually, adding counterterms is equivalent to splitting 
the bare coupling constant into the renormalized coupling constant and some infinite 
additional term. However, using noninvariant regularizations (and, in particular, the 
regularization, breaking the BRST-invariance, which is used here), it is also necessary 
to add counterterms, restoring the Slavnov-Taylor identities [25, 26J in each order of the 
perturbation theory. In general, it is necessary to analyze such counterterms. However, in 
the considered case the situation is simpler. Really, the one-loop two-point Green function 
for the Faddeev-Popov ghosts is finite and does not depend on regularization. Interaction 
of ghosts with the background field is fixed by the background gauge invariance, which 
is unbroken with the considered regularization. Therefore, additional counterterms do 
not contribute to subtraction diagrams, containing a loop of the Faddeev-Popov ghosts, 
in the two-loop approximation. Moreover, terms with the Faddeev-Popov ghosts do not 
evidently depend on whether the bare or renormalized coupling constant is in the gauge 
fixing action. Hence, their contributions do not also depend on a way of splitting the bare 
coupling constant into the renormalized one and counterterms. 

Quantizing the theory we also write the bare coupling constant eo in the gauge fixing 
terms. Therefore, a part of the action, quadratic in the quantum field, is written as 

V. 



^tr Re J d*x d'd V [uj (l + + 2W»D a 



Breaking the invariance under the BRST-transformations can lead to the necessity of 
adding counterterms proportional to 



tr J d A xd 4 6VD^V. (82) 

(If the background field is 0, this follows from Refs. [23, EE]- Terms, containing the 
background field, can be restored from the background gauge invariance.) But this means 



19 



that all one-loop diagrams, including diagrams with insertions of both the counterterms, 
appearing due to the renormalization of the coupling constant, and the additional coun- 
terterms, with a loop of the quantum field V, are 0, because they can contain no more 
than 2 spinor derivatives. 

At last, let us consider diagrams, containing a loop of the Nielsen-Kallosh ghosts. 
Since the Nielsen-Kallosh ghosts exist only in the one-loop approximation, there are 
no additional counterterms, caused by the noninvariance of the regularization under the 
BRST-transformations, in these diagrams. However, the contribution of the counterterm 
diagrams is essential due to the renormalization of the coupling constant. Really, the 
coefficient in the action for the Nielsen-Kallosh ghosts should be the same as in the gauge 
fixing terms. Therefore, it must contain the bare coupling constant: 

-\tr [ d 8 xB + e™B. (83) 
H J 

To regularize diagrams with counterterm insertions and a loop of Nielsen-Kallosh ghosts, 
the action for the corresponding Pauli-Villars fields should be written as 



tr J d'x ( J d'6 B+ v e™B PV + ^§ | dW PV + ^§ | d%B$ v Y) . (84) 

This expression also contains the bare coupling constant e^. Mb is proportional to the 
regularization parameter A. Really, let us present a bare coupling constant as 

\ = \z 3 , (85) 

where e is the renormalized coupling constant, and Z 3 is the renormalization constant. 
Then, expanding the Pauli-Villars determinant for the Nielsen-Kallosh ghosts in powers 
of Z 3 — 1, we obtain terms, regularizing diagrams with insertions of counterterms. 

However, due to inserting this determinant the generating functional starts to depend 
on the normalization point at a fixed bare coupling constant eo, because the renormalized 
coupling constant e depends on ji. 

In the Abelian case calculating divergences for the action, similar to (1541) . was made, 
for example, in Ref. [M]. In the considered case it is also necessary to take into account 
a factor — C2/2, which appears because the Nielsen-Kallosh ghosts are in the adjoint 
representation of a gauge group and anticommute. (There is only one matter superfield 
now, instead of 2 matter superfields in the Abelian case.) Moreover, the renormalization 
constant of a matter field Z should be substituted for the constant Z 3 . Taking into account 
these comments, the result of Ref. [21] can be formulated as follows. Contribution of the 
counterterm diagrams for the Nielsen-Kallosh ghosts to 1/a can be written as 

^lnZ 3 . (86) 

To find this contribution in the two-loop approximation, we note that after the one- 
loop renormalization the renormalization constant will be 

^3 = l + ^m- + 0(a 2 ). (87) 
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Therefore, the contribution of diagrams with counterterm insertions in the two-loop ap- 
proximation is written as 



f^M- (88) 
(2tt) 2 /i 

This contribution exactly cancels the two-loop divergence, so that after the one-loop 
renormalization 

d~\a, pfp) = i- ^ln^-^ln^ + 0(a 2 ). (89) 
a 2ir p {2tc) p 

For an arbitrary order of the perturbation theory it is reasonable to propose that the 
two-point Green function of the gauge field is given by 

1 — -^-C 2 \nd(a ,A/p) + ^-C 2 \nZ 3 (a,A/i2)-^-C 2 \n-. (90) 



d(a,(i/p) a 27r 2ir ' 2ir p 

Really, it is easy to see that the exact NSVZ /5-function is obtained by differentiating 
this equality with respect to Inp, and the term, proportional to lnZ 3 is obtained from 
contributions of diagrams with counterterms insertions. In the two-loop approximation 
this equation agrees with ([79]) . if the contribution of diagrams with counterterm insertions 
is taken into account. 

If Eq. f)90p is true, then divergences exist only in the one-loop approximation. Really, 
because 

1 ' Ma,k/n) (91) 



d(a,p,/p) d(a (a, A///), A/p) 

is finite, it is necessary to cancel only the one-loop divergence. For this purpose the bare 
coupling constant is presented as 

- = - + ftf 2 m-. (92) 

We note that presence of divergences only in the one-loop approximation in this case does 
not mean that the physical /5-function has only the one-loop contribution. Really, the 
physical /3-function is a derivative of the two-point Green function with respect to the 
logarithm of the momentum if proper boundary conditions are imposed. Such function, 
as we already saw, has corrections in all loops. A relation between the divergences and 
the physical /^-function is broken due to the way of the regularization of diagrams with 
the counterterm insertions, which leads to the dependence of the generating functional 
on a normalization point at a fixed bare coupling constant [35]. Thus, similar to the 
electrodynamics, we obtain that the renormalized effective action is exhausted at the 
one-loop, while the Gell-Mann-Low function has corrections in orders of the perturbation 
theory. Note, that this conclusion agrees with Ref. [ID], in which the two-loop /^-function 
for the N = 1 supersymmetric Yang-Mills theory was calculated with the differential 
renormalization |41j . 

So, if Eq. ( |90l) is valid, the Gell-Mann-Low function coincides with the exact NSVZ 
/3-function, and divergences in the two-point Green function exist only in the one-loop 
approximation. 
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4 Schwinger— Dyson equations and Slavnov— Taylor 
identities 



4.1 Schwinger— Dyson equations for the contribution of matter 
superfields 



The calculations, described in the previous section, reveal some interesting features 
of the quantum correction structure, which appear if the higher derivatives are used 
for regularization. In order to partially explain these features, it is possible to use a 
method [121 02] , based on substituting solutions of the Slavnov-Taylor identities into the 
Schwinger-Dyson equations. Here we will discuss only structure of the matter superfields 
contribution to the exact /3-function. Contribution of diagrams with loops of the gauge 
fields and ghosts is not so far calculated using this approach. 

In order to construct the Schwinger-Dyson equations for model ([3]) it is necessary to 
split the action into three parts: the action for the background field, the kinetic term for 
quantum fields, which does not contain the background field, and interaction, in which 
the other terms are included: 



S = S(V) + S 2 (V, 0) + Sj(V, V, 0). (93) 

(Earlier we saw that terms of the first order in the superfield V, which were obtained from 
the expansion of the classical action, can be omitted.) So, generating functional ( 1201) can 
be written as 



Z[J,j, V] = J dfi exp (iS[Y) + iS 2 [V, <p] + iSj[V, V, 0] + iS s + % J d 8 x JV^j = (94) 
= exp (iS[V}+iSj -— x 

dfx exp (iS 2 [V, 4>] + iS s + i J d 8 x JV + i J d 8 x jvj 



where / d 8 x = / d 4 xd 4 9 x . Let us differentiate this expression with respect to the back- 
ground field 



() Z[J,j,V]=^Z : exp (/,S'/ 



5V 



15 15 1 5 1 
i 5 J' i 53 ' i 5j 



zJ x x 



x J dfi exp (iS 2 [V,4>}+iS s + i J d 8 xJV + i J d 8 xJV 



Moving the current J x to the left and dividing the result to Z, we obtain 



5 



S, 



i 5 J i 5 j 



3=0 



(95) 



(96) 



5V X " L " ,J ' " J 5V X ' 5Y X 
Because the background field V is a parameter of the effective action, these equality can 
be equivalently written as 
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5T _ SS[V) i 5 
6VZ ~ 5V X + \ W S 

where the angular brackets denote taking an expectation value by the ordinary functional 
integration. Certainly, it is necessary to set the field V (an argument of the effective 
action) to in the final result. 

Let us find a contribution to this expression given by the matter superfields. The 
corresponding interaction terms are 

Si = — J d 8 x ^e v e 2V 'e v — 1^0 + similar terms with 0. (98) 

Differentiating Si with respect to the background field, we obtain that the corresponding 
contribution to the effective action is written in the form 



(97) 



= Terms without matter + Ci £\Fa ( ^ n ^ y' 1 ^») } 

+-^— j d 8 x (^(j) x ~e Vx e 2Vx e Vx (j) x + similar terms with </y. (99) 

We are interested in the two-point Green function of the gauge field, corresponding to 
the expansion of the effective action in powers of the background field up to the second 
order terms. This function is evidently symmetric with respect to the indexes, numerating 
generators of the gauge group. Then, using the form of action for additional sources fl22|) . 
we easily obtain 



T a e Vx e 2Vx e Vx 



+ 



:e Vx e 2Vx e Vx T a c 



1 

-tr 

i 



rpo 



5 2 r 



5 2 T 



(100) 



(Here the derivatives with respect to the sources must be expressed in terms of fields.) 

Therefore, using Eq. (199]) and taking into account similar terms with the fields </>, the 
corresponding contribution to the two-point Green function of the matter superfield can 
be written as 



5T 



5V15V* 



+ e 



V V y 



tr 



rj-iQ 



5 2 T 



(T 



a\t 



5 2 T 



+ h.c. 



(101) 



1 sj+ s^ox 

where dots denote contributions of the gauge fields, ghosts, and also all possible Pauli- 
Villars fields. We note that the calculation of the Pauli- Villars fields contributions are 
made completely similar to the calculation of ordinary fields contributions [12] , and details 
of this calculation are not presented here. The result will be given below. 

Performing differentiation in Eq. fllOip . this equation can be graphically presented as 
a sum of two effective diagrams 



Ar 



(2) 




(102) 
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The double lines correspond to the effective propagators, which are written as 



5 2 T \ _i GD 2 D 2 X „ B ( 5 2 T \ 1 mJD 



4(«9 2 G 2 + m 2 J 2 ) xyl \<)o.,.<)oJ < p(; - ■ >" 2 -'- XV 



(103) 



depending on the chirality of the ends. The functions G and J are determined by the 
two-point Green functions of the matter superfield as 



PV D 2 D 2 r) 2 F D 2 

dl = ^ G {d 2 )5t „; -^L. = -^ m J(^, (104) 



wherec^ =5 4 (x-y)5 4 (#,-^). 



4.2 Slavnov— Taylor identities and their solutions 



The vertex functions in Eq. (11021) can be obtained from the Slavnov- Taylor identities. 
Certainly, it is necessary to take into account that we consider vertex functions, which 
have an external line of the background field. As we already mentioned above, the 
effective action is invariant under the background gauge transformations. It is easy to see 
that this invariance can be expressed by the equality 



-l-similar terms with and 0o I , (105) 



where 



5<p x = A x (j) x ; 50+ = + A+; 5<f) 0x = A x <p 0x ] 8($ x = (f)^ x A+; 

5V X = ~(A X + At)+0(V). (106) 

Here A is an arbitrary chiral superfield, and all other terms in 5V are proportional at 
least to the first degree of the background field. (For the fields e.t.c. the transformation 
laws can be also easily written.) Let us differentiate Eq. (I105p with respect to Ay, (f>Q Z , 
cj) x or with respect to A" 0q 2 , <p x . As a result we obtain the Slavnov- Taylor identities 



= ^m 2 + J D 2 ) — — — T a D 2 5 8 -D 2 (s 8 T a — 



= ±(D 2 + D 2 ) °-± ±J—T a Dy xy + D 2 (5lT a J^— ), (107) 



where 
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Sr^w- imr^^- 1 ) ^ (108) 



Solutions of Eqs. ( jlOTf) are the functions 



8 3 T 



5V?,t 



J 0z u Vx 



p=0 



2d 2 U 1/2y [Dy xy S 8 yz )F(q 2 ) + \D b C bc Dl{Dl5l y Dl5l)f{q 2 



(109) 



5 3 T 



5V«5(j) 0z 5(j) x 



D 2 D 2 5 8 5 8 ) + 



p=0 



f 2 D b C bc Dl[DlDl5l y Dl5l)h{ q 2 ) + ^V)(%^ 



y y xy yz 



m[J\q 2 ) J(q 2 )-1 



16 



D 2 D 2 

£)2£)2^8 ^2/^2/ £8 , r.2r,2 x 8 „2 X 8 



y y xy 



16 



yz 1 v v xv" • 



y y"xyi "yz 



Here the primes denote derivatives with respect to g 2 

1 



rii/ 2 



D a D 2 D n 



1 D a D 2 D n 



rj-ia 

(110) 

111) 



8<9 2 8<9 2 " 

is a supersymmetric transverse projection operator, and the functions F, f and h can not 
be determined from the Slavnov- Taylor identities. 

Functions (I109p and flllOp allow finding ordinary Green functions by the identities 



d 2 5 3 v s 3 r D 2 6 3 r s 3 r , s 

* — ■ E _ — (112) 

2 5V%5cj )y 5<f>+ SV'StyStf 2 SV'SfaSfa 5V*6<j> y 8$ t 

We note that the additional sources in Eq. (122]) were specially introduced in order that 
such identities take place. Using Eqs. (11121) we find 



5 3 T 



5 3 T 



p=0 

+^GXq 2 )DYl,D y (D 2 y 5l y D 2 y 5l) + ±D 2 y 6 8 xy D 2 y 8 8 yz G(q 2 



d 2 Ui / 2y(D 2 5 8 y D 2 5 8 yz )F(q 2 ) + 



p=0 



- — J'(q 2 ) 
32 W ' 



D 2 5 8 D 2 D 2 5 8 -D 2 D 2 5 8 D 2 5 8 

y^xy y y w yz y y xy y K yz 



rpa 



T a ; (113) 
(114) 



4.3 Exact Gell-Mann— Low function. 



We will calculate 



5T 



dliiA<SVg£Vs 



(115) 



p=0 



25 



Note that the regularization by higher covariant derivatives is essentially used here, be- 
cause it allows differentiating the integrand and taking the limit of zero external momen- 
tum. 

After substituting the vertex functions from Eqs. (jTTSfl , (fill) . ffT09|) . ffTTO]) and the 
propagators from Eqs. ( 11031) . the Weak rotation, and some simple transformations, using 
the algebra of the covariant derivatives, we obtain that in the momentum representation 
the first diagram is written as 



-C(R)- 



d 



-tr 



dlnA 

1 d 
2dq 

and the second one is 



d 4 p d 4 q 



(2tt) 4 (2tt) 4 

+ -(.„( 5 ^ +m ^ 



V9 2 n 1/2 v 

m 2 J 



8GF 



m 2 JJ' 



q 2 G 2 + m 2 J 2 2q 2 (q 2 G 2 + m 2 J 2 ) 
G 2 



q 2 G 2 + m 2 J 2 



q 2 G 2 + m 2 J 2 



+ 

(116) 



-C{R) 



cHnA 



-tr 



d A p d 4 q 



v<9 2 n 1/2 v 



m 2 j(j- r 



(2tt) 4 (2tt) 4 

8Gf + 8m 2 J h 



8GF 



m 2 JJ' 



2q\q 2 G 2 + m 2 J 2 ) q 2 G 2 + m 2 J 2 J q 2 G 2 + m 2 J 2 



V 5 



q 2 G 2 + m 2 J 2 
G 2 



2q 2 (q 2 G 2 + w?J 2 ) 



+ 



;ii7) 



Let do denotes the function d, calculated without taking into account counterterm inser- 
tions on lines of matter superfields, or, equivalently, at Z — 1. Moreover, we take into 
account that for finding the Gell-Mann-Low function it is necessary to set m = 0. Then, 
adding the results for the effective diagrams in the Schwinger-Dyson equation and using 



Eq. (13811 . we obtain 



167rrflnA 



p=0 



MfJp V 



q 2 G 2 PV + M 2 J 2 pv 



+ 



-E c ii( ln (" 2G Pv + A^) + 

i 

^ ( M 2 J PV {2q 2 J> PV -{J PV -l)) 8G PV f PV + 8M 2 J PV h PV \\ 

2q\q 2 G 2 PV + M 2 J 2 PV ) fG 2 PV + M 2 J PV J}' [ * } 

(The dots denote contributions of the gauge fields and ghosts, and the symbol PV means 
that the corresponding function is calculated for the Pauli-Villars fields.) We see that all 
noninvariant terms, proportional to V 2 , and terms, containing the unknown function F, 
are completely cancelled. Nevertheless, there are the functions / and h, which can not be 
found from the Ward identities, in the final result. 

It is important to note that many contributions in this expression are integrals of 
total derivatives. This partially explains the feature, noted earlier. However, calculations 
show that all terms in this expression should be integrals of total derivatives. Moreover, 
the accurate analysis of the calculation, described above, shows that terms, which are 
not factorized into total derivatives in Eq. (1118)) . are always equal to 0. Therefore, it is 
possible to suggest the equality 
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f gg d f m 2 J(2q 2 J'-(J-l)) 8Gf + 8m 2 Jh \_ 
J (27r) 4 dlnAV 2q 4 (q 2 G 2 + m 2 J 2 ) q 2 G 2 + m 2 J 2 J 1 ' 

For the massless theory it can be written in the simpler form 

f d 4 q d f . . 

We will call Eqs. (11191) and (11201) the new identity for the Green functions. It is not a 
consequence of the gauge symmetry, the supersymmetry, or the superconformal symmetry. 
Below we will discuss it in more details, and now let us consider its consequences. If the 
new identity for the Green functions is true, then we obtain 



d 



din A 



d Q 1 {a Q ,A/p) 



16%C(R) 



p=0 



- Ci In (q 2 G 2 PV + M 2 J 2 PV ) - J> 



d A q d Id 
(27r) 4 dlnA V¥ 

M 2 Jp V 
q 2 G 2 PV + M 2 J 2 PV 



\n(q 2 G 2 ) - 



:i2D 



The obtained integral is reduced to the total derivative in the four-dimensional spher- 
ical coordinates, only the substitution at the low limit being different from 



d rfr 1 



din A 



C(R) d 



p=0 2tt din A 



In G(0) 2 - °i ln ( m ?Jpv(0, 2 ) - 



1 1 1 . ./ dlnG 



Jpy(0, Oj) J 7r \ dlnA 



(122) 

g=0 



Here we took into account that the function J P v{q 2 /A 2 , cij) in the limit g — > tended to 
some finite constants, because it was defined by convergent (even in the limit A — > oo), 
dimensionless integrals, which did not contain infrared divergences. Because both parts 
of Eq. (I122p depend on a and A/p, it allows finding the expression for d$ 1 up to an 
insignificant numerical constant 



doVoA/p) = -C(R)( In- -InGW A/p)) + const. (123) 

7T V V J 



A 

p 

Now let us calculate the sum of diagrams with counterterm insertions on lines of matter 
superfields exactly to all orders of the perturbation theory. For this purpose we make the 
substitution 



Q^Q/VZ; <§ -»• $/Vz (124) 

in the generating functional. Then, it is easy to see that if 0o an d j are 0, the dependence 
on the renormalization constant Z can be found by the replacement 

m -> m/Z) Mi -> Mi/Z. (125) 

Making this replacement in Eq. (I122p it is possible to restore the dependence of the 
effective action on Z: 
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d~ 1 (a, n/p) = d,Q 1 (ao, A/p) lnZ(a,A//x) 

7T 



— + -C(fl) f In - - In ZG(a, Li/p)) + const. (126) 
a tc \ p / 



Differentiating this expression with respect to the momentum p, we obtain the Gell-Mann- 
Low function 

amp amp tx V / 

This expression corresponds to a contribution of the matter superfields to the exact NSVZ 
/5-function (pQ). (It is necessary to take into account that in the considered case the matter 
fields are in the representation R + R.) According to Eq. (I126p . to cancel the dependence 
on A, the bare coupling constant should be presented in the form 

— = ---CYi?) In-. (128) 

This means that there are divergences only in the one-loop approximation. 

Note that if the new identity were not valid [13], the contribution of the matter su- 
perfields to the exact /3-function should be modified as follows: 



^w (l _ 7(a) _„ m 

IT \ p^o G 

5 New identity for Green functions 



(129) 



The new identity for the Green functions in the massive and massless cases are given 
by Eqs. (11191) and (11201) respectively. In the both cases it can be equivalently rewritten 
in the following functional form (33] (for the simplicity we will consider here the Abelian 
case) : 

/ d*x d% V y D«V x ; , ; ' 1 =0, (130) 

where we assume that the derivatives with respect to sources should be expressed through 
the derivatives with respect to fields. The condition p = means that in this case the 
background field V depends only on variables 6 and is independent of the usual coordinates 
x^. The derivative with respect to In A in this expression is essential. In order to check 
this, let us consider for simplicity the massless case and suppose that these derivative is 
absent. Then, from the dimensional considerations 

f d A q f r d A q a{q 2 /A 2 ) 
J J2^^G ~ J (2vr)4 ' 

where a is a dimensionless function, which is rapidly decreasing at q 
it is possible that a(0) ^ 0. But if a(0) 7^ 0, then the integral in Eq 



(131) 

* 00. In general, 
ffT3"TT) is not well 
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defined: it is divergent in the infrared region. In order to avoid this we introduce the 
additional differentiation with respect to In A. Due to its presence the term a(0), which 
does not depend on A, disappears, and the integral becomes finite in the infrared region. 
Thus, without the derivative with respect to In A the left hand side of the new identity 
is not well defined. It is essential that the differentiation with respect to In A is possible 
due to using the higher covariant derivative regularization. 

To prove equality (11301) it is necessary to replace differentiation with respect to the 
sources by the differentiation with respect to fields and substitute expressions for the 
vertexes and propagators, presented earlier. 

Let us rewrite identity (11301) in a different form. We commute the differentiations 
with respect to V and 4> , and again use the Schwinger-Dyson equations, keeping only 
contributions of the matter superfields as earlier. The result is 



+0 O V^ - 4>te~ 2Vt 4>) 



= 0. (132) 

z=x,p=0 



A contribution of terms, containing (f)^ is 0, because after performing the differentiation 
we obtain 



d 2 Yx ! d 2 5j; 



= 0. (133) 

lz u V0x z=x 

(In order to verify the last equality it is necessary to express the derivative with respect 
to the source through derivatives with respect to fields, use Eqs. (llOSp and fll03j) . and 

take into account that the expression P X S 4 (8 X — 6 y ) is not only if the operator P 

contains 4 spinor derivatives. 

Therefore, the new identity can be written in terms of composite operators correlators 

as 



d 



din A 



d 8 xD a Vj^cf> + e 2V c 



d»yV y U + e 



,2V A _2+ e -2Vt 



p=0 



■ (134) 



Taking into account the identity 



D 2 D 2 



d 2 



-160 H 



this expression can be rewritten as 



(135) 



d 



din A 



d 8 x D a V T 



D n D 2 



d 2 
D 2 D 2 
—$2— 



e 6 — 



D 2 D 2 ~ 



-2V* 



p=0 



(136) 



Using the Leibnitz rule for the supersymmetric covariant derivative, we find 
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D 2 D 2 D 2 D 2 D D 2 B 

^(V0+) = V^0+ + 2/TV^V - 8i(^C) ab D b D a V^4> + + ..., (137) 

where dots denote all terms, which do not contribute to the considered correlator, because 
the momentum of the field V is 0. (As a consequence, all terms, in which more than 4 
spinor derivatives act on the fields V are 0.) The first term in Eq. (j!37p also does not 
contribute to the considered correlator, because (the sources are not set to 0) 



d\nA\J lQd 2 v ^' V d\nA\J v r ' 16d 2 

d f j8 8 D 2 D 2 5T d f a 5 D 2 5T . . 

rf 8 x r ,, — - = -— — r / rf 8 x — 09r ^, =0- (138) 



din A J 5j+ 49 2 50+ dlnAJ 5j+ 2d 2 
Therefore, the new identity can be equivalently written as 



d 



dlnA 



d 2 



d 8 y 



xD c D b Y„^- 



5 

y ^2 ^2/ 



similar terms with 



0. 



p=0 



(139) 



Actually this equality has been already written in Ref. [44] on the language of Feynman 
diagrams. Nevertheless, its strict functional formulation is first presented here. 



6 Verification of new identity for Green functions 
6.1 Four- loop approximation in the Abelian case 



In order to be sure that the calculations, presented above, and the proposal about 
existence of the new identity for Green functions are true, it is necessary to perform a 
verification by explicit calculations. Making this verification is considerably simplified, if 
we note that expressions (11161) and (11171) allow not only finding sums of all diagrams with 
two external lines of the gauge field, but also summing special classes of such diagrams. 
Such classes of diagrams are obtained from a frame, to which external lines are attached 
by all possible ways. In particular, in order to verify Eqs. (11 16ft and (11171) in Ref. [15] we 
considered a group of diagrams, which were obtained from a frame, presented in Fig. [21 
by attaching two external line of the gauge field. 



Figure 2: Schematic picture of group of diagrams, which is used for verification of new 
identity 
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This diagrams can be calculated by two different ways: 

1. using Eqs. (11161) and (11171) and the functions G, f and F, obtained preliminary 

2. by explicit calculation using supergraphs. 

Thus it is possible to perform a four-loop verification of Eqs. (11161) and (I117p . and 
also identity ffT20l) . 

In order to find diagrams which will be essential for obtaining the unknown functions 
G, f and F in the considered case, it is convenient to use the following simple speculations: 
A diagram, presented in Fig. [21 can be considered as a formal product of the one- and 
two- loop diagrams with <p* and <fi external lines (pairs of their ends are identified): 



/0\ 




or as the three-loop diagrams with the identified ends: 




The parts of the diagram, obtained by this way, are the Feynman diagrams for finding 
the function G. To find the functions / and F, it is necessary to add one more external 
line of the superfield V to such diagrams. 




Figure 3: The diagrams, used for calculation of the function G 

Thus, it order to find the function G in the considered case it is necessary to calculate 
diagrams presented in Fig. [21 Then the function G is obtained according to definition 
( 11041) . One-, two- and three-loop parts of the function G, defined by diagrams in Fig. 
El will be denoted by Gi, G2 and G3 respectively. Expressions for them, obtained by 
calculating these diagrams, are 

, k 2n \ -1 1 

p^( 1+ a^) mrw ; (140) 
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G2= _4eM^fi + sir'fi + ' 2 " rI {2q+k+l) 



(2tt) 8 V k 2n ) V k 2n ) k 2 P(k + q) 2 (q + l) 2 (k + q + l)^ 

(141) 

« , d 4 kdHd A r ( P n \-i/ / 2n \-i/ r 2 " x -i 
G 3 = -4e b — — — 1 + -^- h+_) h + 



+ 



(2tt) 12 V A 2 "/ V A 2 "/ V A 2n , 
1 / (2q + 2k + r + lf{q + I) 2 



k 2 l 2 r 2 (q + k) 2 (q + l) 2 (q + k + I) 2 \(q + r + k) 2 (q + r + k + I) 2 

(2q + k + l) 2 2{2q + k + l) 2s 
(q + r + k + l) 2 (q + r + kf 

The complete function G (certainly without diagrams, which are not essential for this 
paper) in the considered approximation is given by 



G(q 2 ) = 1 + G 1 + G 2 + G 3 , (143) 

where the unity is a tree contribution. 

The functions F(q 2 ) and f{q 2 ) are defined from three-point Green functions by Eq. 
( 11091) . As in the case of the function G, in order to obtain a contribution, which corre- 
sponds to the considered class of the four-loop diagrams, it is sufficient to calculate only 
diagrams of a special form. They are obtained from diagrams, presented in Fig. [3X by all 
possible insertions of one external ^/-line. In all these diagrams one external straight line 
corresponds to the chiral field <fi and the other corresponds to the non-chiral field 0q. 

We will denote one-, two- and three-loop contributions to the function / by fi, f 2 
and f 3 respectively. Similar notation we will use for the function F. Because in the tree 
approximation these functions are and f\ — 0, in the considered order we have 

F = F 1 + F 2 + F 3] f = f 2 + f 3 . (144) 

Expressions for fx, f 2 , f 3 , F\ and F 2 , obtained by calculation of the above pointed diagrams 
have the following form: 

A(g 2 ) = 0; (145) 

1 f d 4 k d 4 l / l 2n 



2 J (2vr) 8 V A 2 V V A 2 " 
1 

X k 2 I 2 (k + q) 2 (q + l) 2 (k + q + I) 2 * 

/ (2q + k + l),(k + qy (2q + k + l),{l + qy \ 

V ik + q) 2 + (l + q) 2 Y {1Ab) 

C d 4 k d 4 l d 4 r / k 2n \ ~ l / l 2n \ - 1 / r 2n \ - 1 
A(9 . ) = 4e -/^I (l + ^) (l + ^) x 

1 f I 1 



k 2 I 2 r 2 (k + q) 2 (q + I) 2 (q + k + I) 2 \(k + q + r) 2 2(k + I + q + r) 2 
(2q + k + iy(k + qY {2q + k + iy{l + qy {2q + k + iy{k + qy 



2(k + q) 2 (k + q + r) 2 2(1 + q) 2 (k + q + r) 2 2(k + q) 2 (k + I + q + r) 2 
{2q + k + iy{k + q + ry s 
2{k + q + r) 4 



(147) 
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e 2 f d 4 k ^ k 2n \-^ 



F ^> = ~jJW)->[ 1 + w«) ww^W 1 (148) 



4 J (2vr) 8 V AW V A 2n 

^4-2#±£-^ + ,* + 5). (149) 



k 2 l 2 (k + q) 2 (q + l) 2 (k + q + l) 2 \ (k + q) 2 (q + k + I) 

An expression for F 3 has been calculated, but it is not presented because it is very large. 
(It is not required for verification of identity ( 11201) .) 

Using the obtained expressions it is possible to verify identity (I120p . With the consid- 
ered accuracy we have 



./ 



(/s - /a • Gi) . (150) 



q 2 G q 2 

Substituting here expressions for the functions G\, /2 and / 3 from Eqs. ()140p . (I146p . 
(11471) we obtain, that the integrand can be written as a total derivative with respect to 
the momentum q 



d A q d f 6 f d 4 q d A k dH d A r d 



(2tt) 4 din A q 2 G J (2tt) 16 din A 

2n 

X( 1 + ' 



1 + X 



\2n J \ J{2n 



/-'■N--' d f {2q + k + iy 



A 2r, 



dq» { k 2 l 2 r 2 q 2 (k + q) 2 (q + l) 2 (q + r) 2 (k + q + I) 2 



(151) 



This equality can be checked by calculating the derivative with respect to q^ using the 
Leibnitz rule and comparing the result with Eq. (I150p . in which expressions for the 
functions G\, j'2 and / 3 are obtained by explicit calculation of diagrams. 

Because the integrand in Eq. (I15ip is a total derivative with respect to q^ of the 
expression, which goes to in the limit q — > 00, expression (11511) is 0. This means that 
identity (11201) is correct in the considered approximation and for the considered class of 
diagrams. 

Note that the equalities, present above, allow checking the method of summing Feyn- 
man diagrams by using the Schwinger-Dyson equations and Ward identities. For this 
purpose it is possible to calculate both effective diagrams in Eq. (11021) in the four-loop 
approximation explicitly and compare the result with Eqs. f II 1 6 f) and (I117p . In considered 
approximation 

?5? k(G2) = 7^( G *- G '- G *> (152) 

\(F 3 -F 2 -G 1 -F 1 -G 2 y (153) 



q 2 G q 

Explicit calculations were made by the method, proposed in Ref. [16], which simplified 
finding a part of diagram, proportional to Vd 2 Ui/2V. Nevertheless, it does not allow to 
obtain a part, proportional to V 2 . That is why the verification of Eqs. (I116p and (11171) 
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was made only for terms, proportional to Vd 2 Ui/2V. In the both cases this verification 
completely confirms them. 

As a small technical remark let us note, that making this verification it is necessary 
to take into account that a large number of ordinary Feynman diagrams contributes both 
to the first effective diagram in Eq. (11021) and to the second one. For example, it is easy 
to see, that a contribution of a diagram, presented in Fig. 0] 




Figure 4: One of the diagrams, giving contribution to both effective diagrams 

is divided into parts 3/4 and 1/4, which correspond to the first and to the second diagrams 
in Eq. ffT02l) . 

6.2 Three-loop approximation in the non-Abelian case 

The Feynman rules are different in a non-Abelian theory mostly due to vertexes with 
the selfaction of the gauge field. That is why it is desirable to check the new identity 
also in this case [47J. For diagrams that do not contain such vertexes the calculations are 
similar to the Abelian case. But for diagrams with the triple vertex of the gauge field 
an additional verification is very desirable. As we already mentioned, for this purpose it 
is not necessary to calculate all Feynman diagrams in a given order of the perturbation 
theory. It is sufficient to consider, for example, a typical three-loop diagram, presented 
in Fig. [5] 




Figure 5: Diagram, generating considered contribution to the new identity 




Figure 6: Way of cutting the diagram 

From the topological point of view there is the only way to cut a loop of the matter 
superfield, presented in Fig. O Hence, it is necessary to calculate a set of diagrams, 
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presented in Fig. [7J As earlier, in all these diagrams the chiral field <fi is at the first 
external line and the non-chiral field 0q is at the second line. Therefore, all presented 
diagrams are not topologically equivalent. 




Figure 7: Diagrams, defining the function /, corresponding to the considered diagrams 

Calculating these diagrams we can find the function /. The function G in the lowest 
approximation should be set to 1. Really, in the tree approximation G = 1. Hence, in the 
given order for the considered class of diagrams we have: 

G{q 2 ) = 1 + 0(a 2 ); f(q 2 ) = f 2 (q 2 ) + 0(a 3 ). (154) 
where f 2 is proportional to a 2 . Therefore, 



J4„ d f(<l 2 ) [ A„ d fM , r.fR-s 



d\nAq 2 G(q 2 ) J din A q 2 

So, we see that the considered contribution is actually determined by the two-loop value 
of the single function / 2 . 

In order to find the function f 2 in the two-loop approximation, it is necessary to make 
an explicit calculation of Feynman diagrams, presented in Fig. [3, using the standard 
supergraph technique. The result is (in the Euclidean space, after the Weak rotation) 

/,(,) = -Wft^C*) - \C 2 ) j ^(jj-^, + x (156) 

x (k + q + l)^ 

(k + q) 2 {k + q + l) 2 k 2 (l + k 2n /A 2n ^j I 2 (l + l 2n / A 2n ^j (k + I) 2 (l + (k + l) 2n /A 2n ^j ' 

where C 2 (R) and C 2 are defined by 
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T a T a = C 2 {R), 



amn fbmn 



jamn j 



Co 5' 



ab 



(157) 
(158) 



Substituting this expression into Eq. (I155p . we obtain that in the considered approxima- 
tion for the considered diagrams 



A d f(q 2 



(2vr) 4 d\nAq 2 G(q 2 ) 
a 2 TT 2 C 2 (c 2 (R)-~C 2 ) J 



d q a Ik d I d \ , d 
A— 



(2vr) 12 dqv { dA 



(k + q + iy 



x 



x 



k 2 (l + fc 2n /A 2n ) I 2 (l + / 2 "/A 2n ) (k + I) 2 (l + (k + l) 2n /A 2n ^j 



q 2 (k + q) 2 (k + q + l) 2 

= 0. (159) 



This means that the new identity for Green functions is also valid in the non-Abelian 
theory. 



7 Conclusion 



Our investigation of applying the higher derivative regularization to calculation of 
quantum corrections in supersymmetric theories allows revealing some interesting fea- 
tures, which were not noted earlier. We should first mention the completely unexpected 
result: all complicated integrals, defining the Gell-Mann-Low function, are integrals of 
total derivatives and can be easily calculated. It is important to note that with the higher 
derivative regularization the Gell-Mann-Low function is calculated in the simplest way. 
According to our calculations this is a function that coincides with the exact NSVZ (3- 
function. So, it is not necessary to tune a renormalization scheme, as it was done earlier, 
calculating the /^-function defined by divergence in the MS'-scheme with the dimensional 
reduction. Note that with the higher covariant derivative regularization there are diver- 
gences only in the one-loop approximation. In a certain degree this confirms speculations, 
presented in Ref. [37|, in which the authors suggested that the Wilsonian action was 
exhausted at the one-loop. The calculations show that with the higher derivative regu- 
larization the renormalized action is exhausted at the one- loop. 

The features, which appear calculating quantum corrections in supersymmetric theo- 
ries, can be partially explained by substituting solutions of the Slavnov-Taylor identities 
into the Schwinger-Dyson equations. However, it is necessary to suppose existing a new 
identity for the Green functions, which does not follow from known symmetries of the 
theory. The exact /5-function can be also considered as a similar identity. Quite possible 
that all such identities are consequences of some nontrivial symmetry. Existence of such 
symmetries was already proposed earlier [TT] in iV = 1 finite supersymmetric theories. 
That is why the obtained results could be also useful for their investigation. In particu- 
lar, it is possible that the anomalous dimension in iV = 1 finite theories is also an integral 
of a total derivative. Now this statement is being checked by the explicit calculation. 
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Finally, it is necessary to enumerate some interesting problems, which have not yet 
been solved. One of them is a proof of the new identity starting from some symmetry. May 
be, there is a relation between the considered problem and the AdS/CFT-correspondence. 
In favour of this we note that it is convenient to formulate the new identity in terms of 
vacuum expectation values of gauge invariant operators. Possibly it would be possible to 
derive the NSVZ /5-function for the pure Yang-Mills theory exactly to all orders of the 
perturbation theory. In this case the Schwinger-Dyson equations are rather involved and 
their investigation is much more complicated. 

To conclude, we can say that dynamics of supersymmetric theories is unexpectedly 
rich by the interesting features, which sometimes are hard to be explained. The higher 
covariant derivative regularization is an excellent tool for revealing these features. 
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